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FOREWORD 


This  report  along  with  the  previously  de¬ 
livered  data  diskette  and  Technical  Information  Memo 
randura  (Ref.  1)  have  been  prepared  in  response  to 
needs  identified  at  the  Gravity  Gradiometer  Survey 
System  (GGSS)  Data  Analysis  Working  Group  meeting  of 
4  April  1984.  The  work  described  herein  was  per¬ 
formed  for  the  Air  Force  Geophysics  Laboratory  and 
the  Defense  Mapping  Agency  under  Contract  No. 
F19628-83-C-0053 ,  "Independent  Analysis  of  Gravity 
Gradiometer  Survey  System  Development  and  Deployment 
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1. 


INTRODUCTION 


1.1  OVERVIEW 

This  report  addresses  the  problem  of  processing  air¬ 
borne  gravity  gradiometer  measurements  to  estimate  the  gravity 
disturbance  vector  at  the  surface  of  the  earth.  The  approach 
is  based  on  a  minimum-variance  estimation  technique,  which  is 
optimal  with  respect  to  probabalistic  models  for  the  uncer¬ 
tainties  associated  with  the  measurement  errors  and  the  a  priori 
statistical  information  about  the  gravity  disturbance  vector. 

A  primary  result  of  the  work  reported  herein  is  a 
practical  methodology  for  processing  the  gradiometer  measure¬ 
ments  to  meet  specified  accuracy  requirements.  The  methodol¬ 
ogy  consists  of  two  stages.  The  first  stage  involves  averaging 
the  individual  gradiometer  measurements  over  appropriate  sub- 
regions  within  the  survey  area  using  a  data  template.  During 
the  second  stage,  these  averaged  gradiometer  measurements  are 
optimally  weighted  and  summed  to  estimate  the  three  components 
of  the  gravity  disturbance  vector  at  the  surface  of  the  earth. 
The  weights  are  chosen  optimally  with  respect  to  probabilistic 
models  for  (1)  the  measurement  errors  and  (2)  the  a  priori 
uncertainty  about  the  true  gravity  disturbance  vector.  The 
rms  accuracy  of  the  estimated  disturbance  vector  is  character¬ 
ized  by  an  error  covariance  matrix. 

In  preparation  for  testing  of  the  Gravity  Gradiometer 
Survey  System  (GGSS),  a  prototype  example  of  an  averaging 
template  is  documented  in  this  report.  The  gravity  and  meas¬ 
urement  error  covariance  matrices  (valid  for  this  template  and 
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the  Attenuated  White-Noise  (AWN)  (Ref.  2)  worldwide  gravity 
model)  are  provided  in  Ref.  1.  This  template  provides  a  simple 
but  modestly  accurate  example  of  the  measurement  averaging 
portion  of  the  data  processing  methodology.  Other  templates, 
having  more  averaging  zones  for  meeting  the  required  rms  accu¬ 
racy  and  optimized  for  data  processing,  are  being  developed  at 
TASC  specifically  for  GGSS  testing.  The  prototype  template 
described  here  is  provided  as  an  interim  example  of  the  type 
of  processing  that  will  be  required. 


1.2  TECHNICAL  APPROACH 

The  technical  approach  of  this  study  is  based  on 
minimum-variance  estimation  theory.  This  approach  has  the 
flexibility  needed  to  treat  mixtures  of  different  types  of 
gravity  quantities  acquired  at  different  survey  heights.  The 
approach  applies  to  both  point  and  mean  gravity  field  data. 

The  central  aspect  of  this  study  is  the  design  of  the 
data  averaging  templates.  The  templates  are  being  developed 
to  achieve  required  rms  accuracy  levels.  This  is  accomplished 
by  using  error  covariance  analyses  to  study  the  tradeoffs  be¬ 
tween  estimation  accuracy  and  the  template  parameters  for  dif¬ 
ferent  blends  of  gradient  data.  The  covariance  analyses  can 
be  performed  before  any  survey  data  are  available  because  the 
error  covariances  are  computed  using  gravity  and  measurement 
e  r ro  r  mode  1 s . 


1.3  ORGANIZATION  OF  REPORT 


The  theoretical  basis  for  the  estimation  algorithm  is 
presented  in  Chapter  2.  First,  the  quantities  to  be  estimated 


L 


and  the  quantities  to  be  measured  are  defined.  The  practical 
need  for  the  data  processing  to  reduce  the  dimensionality  of 
the  computations  is  explained.  Based  on  these  preliminaries, 
the  estimation  problem  is  stated  and  a  recursive  algorithm  for 
solving  it  is  presented. 

The  remainder  of  Chapter  2  describes  some  recommended 
conventions  for  documenting  the  data  analysis.  A  prototype 
data  processing  template  is  defined  and  the  derivation  of  the 
covariances  which  are  needed  to  design  an  optimal  estimator 
for  this  template  is  explained.  The  report  ends  with  Chap¬ 
ter  3 ,  which  presents  a  summary  and  conclusions. 
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2. 


TECHNICAL  DISCUSSION 


The  purpose  of  this  chapter  is  to  formulate  and  dis¬ 
cuss  the  problem  of  processing  data  from  an  airborne  gravity 
gradiometer  survey.  The  principal  result  of  this  analysis  is 
a  practical  algorithm  for  (1)  estimating  the  point  gravity 
disturbance  vector  at  the  surface  of  the  earth  (or  other  appro 
priately  defined  downward  continuation  point)  from  the  survey 
data  and  (2)  computing  mean-square  errors  of  those  estimates. 

2.1  STATEMENT  OF  THE  PROBLEM 

Processing  survey  data  to  estimate  the  gravity  dis¬ 
turbance  vector  is  formulated  as  a  minimum-variance  estimation 
problem.  This  formulation  involves  three  important  elements: 

•  Definitions  of  the  gravity  quantities 
that  are  being  estimated 

•  Definitions  of  the  gravity  quantities 
that  are  being  measured 

•  Probabilistic  models  for  the  uncertainties 
associated  with  the  survey  measurement 
data  and  all  associated  gravity  quanti¬ 
ties  (both  those  being  measured  and  those 
being  estimated). 

2.1.1  Estimated  Gravity  Quantities 

The  gravity  quantities  to  be  estimated  depend  on  the 
purpose  of  the  data  processing.  For  data  processing  that  will 
follow  routine  GGSS  survey  operations,  the  gravity  disturbance 
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vector  at  the  surface  of  the  earth  will  typically  be  estimated. 
In  contrast,  for  the  verification  tests  of  the  GGSS ,  a  spe¬ 
cially  defined  gravity  disturbance  vector  should  be  estimated, 
as  explained  in  the  following  discussion. 

The  gravity  disturbance  vector  at  position  r  is  de¬ 
fined  as  follows: 


6(r)  4  g(r)  -  gR(r)  (2.1-1) 

In  Eq .  2.1-1,  g(r)  is  the  whole  value  gravity  vector  at  r,  and 
gR(r)  is  the  reference  gravity  vector.  The  reference  field 
and  the  coordinate  systems  for  expressing  gravity  vectors  and 
position  vectors  must  be  specified  to  complete  the  definition 
of  the  disturbance  vector.  Appropriate  candidate  reference 
systems  are  WGS  72  or  WGS  84. 

Although  6(r)  is  the  quantity  that  is  usually  sought, 
it  is  inappropriate  to  estimate  this  gravity  disturbance  vector 
during  airborne  testing  of  the  GGSS.  The  reason  is  that  ini¬ 
tial  gradiometer  survey  testing  will  be  confined  to  a  limited 
geographic  area  (probably  on  the  order  of  500  km  by  500  km). 
This  limited  extent  makes  long-wave length  testing  of  the  survey 
accuracy  imprecise.  However,  the  medium-  and  short -wavel ength 
accuracy  of  survey  data  can  be  tested  by  estimating  a  redefined 
gravity  disturbance  vector,  which  in  this  report  is  called  the 
residual  gravity  d  i  s  tyi  rbance  .  This  residual  vector  is  defined 
so  that  its  variance  at  long  wavelengths  can  be  controlled 
to  suit  the  accuracy  objectives  of  GGSS  testing.  The  residual 
gravity  disturbance  vector,  denoted  d(r),  is  defined  as  the 

departure  of  the  disturbance  6(r)  from  its  local  mean  6  (r): 

— m  - 

d ( r  )  A  6 ( r  )  -  6  ( r  )  (2.1-2) 
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For  this  style  of  template,  which  is  comprised  of  con 

t  h 

centric  circuits,  the  documentation  matrix  D  is  n*8.  The  kL 
gravity  quantity,  whose  value  is  given  by  x(k) ,  is  described 
by  the  following  eight  numbers  in  the  kt^1  row  of  D: 

t  h 

D(k,l)  =  ID  Number  of  the  kL  gravity  quantity 

D(k,2)  =  x  coordinate  (xc)  of  the  upper-left  corner 
of  the  template  circuit  (km) 

D(k,3)  =  y  coordinate  (y  )  of  the  upper-left  corner 
of  the  template  circuit  (km) 

D(k,4)  =  width  (w)  of  the  rectangular  subregions  (km) 

D(k,5)  =  mean  altitude  (A,)  of  rectangular  subregion 
No.  1  (km) 

D(k,6)  =  mean  altitude  (A~)  of  rectangular  subregion 
No.  2  (km)  Z 

D(k,7)  =  mean  altitude  (A9)  of  rectangular  subregion 
No.  3  (km)  J 

D(k,8)  =  mean  altitude  (A,)  of  rectangular  subregion 
No.  A  (km) 

For  a  point  gravity  quantity  at  the  template's  center,  as  dis¬ 
tinguished  from  a  mean  quantity,  the  template  circuit  reduces 
to  a  point  at  the  origin,  which  yields  D(k,2)  =  D(k,3)  = 
D(k,A)  =  0  and  D(k,5)  =  D(k,6)  =  D(k,7)  =  D(k,8).  For  the 
mean  z  gravity  disturbance  (i.e.,  when  ID  No.  =  2),  the  mean 
refers  to  an  average  over  a  single  square  area  centered  at  the 
origin.  Therefore,  if  ID  =  2,  the  parameter  w  denotes  the 
length  of  each  side  of  this  square  area. 

Other  forms  of  templates,  having  different  symmetries 
are  appropriate  for  estimating  x  (east)  and  y  (north)  gravity 
disturbance  components  (i.e.,  the  deflections  of  the  vertical) 
Such  templates  are  being  developed  and  will  be  provided  in 
subsequent  reports. 


discussion,  each  of  the  picture- frame  averaging  zones  will  be 
referred  to  as  a  "circuit."  For  mathematical  tractabi 1 i ty , 
each  of  these  circuits  is  composed  of  simple  rectangular  sub- 
regions.  The  recommended  way  of  defining  a  single  circuit 
using  four  rectangular  subregions  is  illustrated  in  Fig.  2.4-2. 
The  geometry  of  each  circuit  is  then  completely  specified  by 
seven  parameters:  the  coordinates  (xc>yc)  of  the  upper-left 
corner,  the  width  w  of  each  rectangular  subregion,  and  the  alti¬ 
tude  of  each  rectangular  subregion.  The  width  w  represents  the 
crosstrack  resolution  of  the  gradient  measurements,  given  that 
the  actual  aircraft  paths  are  expected  to  have  minor  deviations 
from  straight  paths.  a  11327 


SUBREGION 
NO.  1 


A  Typical  Template  Circuit,  Specified  by  the 
Three  Parameters  xc ,  y  ,  w  and  Altitude  of 

Each  Subregion  (Not  to  Scale) 


Figure  2.4-2 


As  discussed  in  Section  2.1,  the  survey  data  are  proc¬ 
essed  by  averaging  them  over  subregions  of  the  surveyed  area. 

(It  is  expected  that  the  data  processing  will  include  track¬ 
crossing  adjustments  to  reduce  low-frequency  measurement  errors.) 
The  resulting  mean  gradient  data  are  listed  in  the  gravity 
measurement  vector  z^  .  According  to  Eq .  2.1-7,  these  averaged 
data  are  a  linear  combination  of  components  of  the  truth  vec¬ 
tor  x  and  the  measurement  error  vector  v^  .  To  document  the 
gravi ty- re lated  quantities  in  x,  an  n*p  documentation  matrix  D 
is  defined,  where  n  is  the  number  of  quantities  listed  in  x, 
and  p  is  the  number  of  parameters  used  to  describe  each  of 
these  quantities.  The  parameters  depend  on  the  particular 
template  that  is  used  for  the  data  averaging,  as  discussed  in 
the  following  section. 

2 . 4  DATA  TEMPLATES 

For  concreteness,  the  discussion  of  data  templates 
focuses  on  a  simplified  example  of  practical  importance  in  test¬ 
ing  the  GGSS:  estimating  the  z  (up)  component  of  the  residual 
disturbance  vector  using  gradiometer  measurements  of  the  zz 
element  of  the  gravity  gradient.  An  appropriate  template  for 
averaging  these  gradiometer  measurements  would  have  azimuthal 
symmetry.  This  symmetry  applies  to  data  processing  that  is 
optimal  with  respect  to  gravity  models  having  an  isotropic 
crosscorrelation  function  relating  the  vertical  disturbance 
with  the  measurements  of  the  zz  gradient  element.  Thus,  a 
reasonable  template  with  complete  azimuthal  symmetry  would 
have  concentric  annular  averaging  zones  like  a  bulls-eye  tar¬ 
get.  However,  airborne  gradiometer  surveys  will  yield  gradient 
measurements  taken  along  nearly  straight  tracks  that  form  a 
rectangular  grid.  To  be  consistent  with  this  survey  geometry, 
the  template  should  also  have  a  rectangular  geometry.  There¬ 
fore,  a  reasonable  template  geometry  for  reduction  of  test 
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2.3 


GRAVITY  DATA  DOCUMENTATION 


The  analysis  of  gravity  gradiometer  survey  data  in¬ 
volves  several  kinds  of  gravity  quantities.  These  include 
point  values  of  gravity  disturbance  vectors  and  disturbance 
gradients,  expressed  with  respect  to  a  reference  field.  Other 
important  gravity  disturbance  quantities  are  mean  values  of 
gravity  disturbance  vectors  and  disturbance  gradients,  averaged 
over  template  zones.  For  convenience  in  designing  future  sur¬ 
veys  and  for  storing,  processing,  and  documenting  actual  survey 
data,  each  gravity  quantity  can  be  assigned  an  identification 
(ID)  number.  These  numbers  are  selected  sequentially.  The 
beginning  of  such  a  list  follows;  additional  quantities  can  be 
added  to  the  list  as  appropriate. 


ID  NUMBER  GRAVITY  QUANTITY  UNITS 


1 

Point 

z 

(Up)  Disturbance 

mgal 

2 

Mean 

z 

(Up)  Disturbance 

mgal 

3 

Point 

X 

(East)  Disturbance 

mgal 

4 

Mean 

X 

(East)  Disturbance 

mgal 

5 

Point 

y 

(North)  Disturbance 

mgal 

6 

Mean 

y 

(North)  Disturbance 

mgal 

7 

Point 

zz 

Disturbance 

Gradient 

mgal/km 

8 

Mean 

zz 

Disturbance 

Gradient 

mgal/km 

9 

Point 

XX 

Disturbance 

Gradient 

mgal/km 

10 

Mean 

XX 

Di s  turbance 

Gradient 

mgal/km 

11 

Point 

xz 

Disturbance 

Gradient 

mgal/km 

12 

Mean 

xz 

Disturbance 

Gradient 

mgal/km 

13 

Point 

yy 

Disturbance 

Gradient 

mgal/km 

14 

Mean 

yy 

Disturbance 

Gradient 

mgal/km 

15 

Point 

xy 

Disturbance 

Gradient 

mgal/km 

16 

Mean 

xy 

Disturbance 

Gradient 

mgal/km 

17 

Point 

yz 

Disturbance 

Gradient 

mgal/km 

18 

Mean 

yz 

Disturbance 

Gradient 

mgal/km 

(entries  to  be  added  as 
future  needs  arise). 
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2.2 


ESTIMATION  ALGORITHM 


The  solution  to  the  estimation  problem  of  the  previous 
section  is  given  by  the  following  recursive  equations  (Ref.  3), 
which  update  the  prior  estimates  and  their  error  covariances. 
The  optimal  estimate  of  the  truth  vector  x  and  the  error  covari¬ 
ance  matrix  for  this  estimate  are 


-1  =  -0  +  K1(z1  -  H-^Xq  ) 


(2.2-1) 


P!  =  P0  ‘  K1H1P0 


(2.2-2) 


In  Eqs .  2.2-1  and  2.2-2,  the  nxm  gain  matrix  is 
computed  as  follows: 


K1  =  P0H!|H1P0H1  +  ri1_1 


(2.2-3) 


The  estimate  of  the  residual  disturbance  vector  and 

the  covariance  matrix  P,.  of  the  errors  in  this  estimate  are 

dd 

given  by  the  following  two  formulas: 


d(r)  =  Bx^ 


(2.2-4) 


Pdd  =  BP1B 


(2.2-5) 


For  the  case  where  the  prior  estimate  Xq  =  0,  the  optimal  esti¬ 
mate  can  be  written  directly  in  terms  of  the  averaged  data  z^ 
as  follows: 


d(r)  =  Gz^ 


(2.2-6) 


G  =  BK, 


(2.2-7) 
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2.1. A  Estimation  Problem 


The  problem  of  estimating  the  residual  gravity  dis¬ 
turbance  vector  d(r)  from  averaged  gradiometer  measurements 
can  be  stated  as  follows: 


Given :  1.  The  a  priori  estimate  Xq  of  the  truth 

vector  x. 

2.  The  covariance  matrix  Pq  of  the  errors  in 
this  prior  estimate. 

3.  The  measurement  vector  of  averaged  gradi 
ometer  data. 


A.  The  covariance  matrix  R-^  of  the  errors  in 
these  averaged  data. 

Find:  1.  An  estimate  d(r)  of  the  residual  gravity  dis¬ 

turbance  vector  d(r)  at  the  template's  center. 


2. 


The  covariance  matrix  P  ,  ,  of  the  errors 

dd 

3(r)  =  d(r)  -  d(r) . 


Optimality :  The  estimates  are  to  be  optimal  in  the  sense 

that  the  weighted  mean-square  errors,  J  and 
,  are  minimized  for  any  positive-definite 
symmetric  weighting  matrices  and  : 


J 

x 


4 


w 


x 


(2.1-20) 


Jd  4  E(3(r)T  Wd  3(r) ] 


(2.1-21) 
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B  =  [1  -1  0lxl2] 


(2.1-12) 


H1  =  l012x2  112xl21 


(2.1-13) 


An  important  conclusion  to  be  drawn  from  this  discussion  is  that 
the  definitions  of  the  truth  vector  x,  the  observation  matrix 
,  and  the  output  matrix  B  are  each  selected  to  suit  the  needs 
of  the  particular  estimation  problem  that  is  being  formulated. 


The  prior  knowledge  of  the  truth  vector  x  is  modeled 
by  an  unbiased  estimate  denoted  Xq : 


x0  =  E(x  ] 


(2.1-14) 


Usually  Xq  =  0  because  a  reference  field  has  been  subtracted  from 
all  gravity-related  quantities.  The  uncertainty  of  this  prior 
estimate  is  modeled  by  its  error  covariance.  The  error  vector 
Xq  and  its  nxn  covariance  matrix  Pq  are  defined  as  follows: 


*0  $  50  '  ^ 


(2.1-15) 


P0  £ 


(2.1-16) 


The  prior  knowledge  of  the  measurement  error  vector 
Vj  is  modeled  by  stating  that  its  mean  is  zero,  its  covariance 
matrix  is  (a  known  matrix),  and  its  crosscovariance  with 


the  error  vector  Xq  is  zero: 


Elvx)  =  0 


(2.1-17) 


EfVjvJ]  =  Rx 


(2.1-18) 


E ( Y1Xq ]  =  0 


(2.1-19) 
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rf.r 


The  truth  vector  in  Eq.  2.1-9  serves  only  as  an  illustrative 
example  of  how  x  can  be  organized;  this  particular  blend  of 
gravity  quantities  would  not  actually  be  used  in  practice  be¬ 
cause  additional  gradient  quantities  (besides  T  )  must  be 
included  to  estimate  the  residual  gravity  disturbance  vector. 

For  the  example  case  under  discussion,  the  number  of 
measurements  in  z^  is  m  =  12  and  the  number  of  gravity-related 
scalars  in  x  is  n  =  18.  The  B  matrix  in  Eq.  2.1-9  is 


B  ~  ll3*3  ~13*3  °3*12  ^ 


(2.1-10) 


where  13x3  is 
zero  matrix. 


the  3x3  identity  matrix,  and  ^3*12  t^ie  ^*12 
This  zero  matrix  accounts  for  the  fact  that  the 


12  gradient  quantities  do  not  enter  the  definition  of  the  point 


residual  disturbance. 


To  simplify  the  bookkeeping,  the  values  of  the  mean 
gradients  in  x  are  listed  in  the  same  order  as  their  measured 
values  in  z  ^ .  With  this  convention,  the  observation  matrix 
has  the  following  form: 


H1  =  [012x6  I12xl21 


(2.1-11) 


In  Eq.  2.1-11,  the  12x6  zero  matrix  indicates  (1)  that  there 
are  12  scalar  quantities  in  the  measurement  vector  z^ ,  and  (2) 
that  the  first  six  elements  of  x  are  not  being  measured. 


Later  in  this  report,  a  prototype  template  is  defined 
for  estimating  the  z  (up)  component  of  the  residual  gravity  dis¬ 
turbance  dz(r),  rather  than  all  three  components  of  d(r). 

For  this  simplified  case,  the  truth  vector  x  is  14x1.  The 
first  two  elements  of  x  are  the  scalars,  6z(r)  and  its  local 
mean,  and  the  B  and  matrices  are 
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d(r) 


Bx 


(2.1-8) 


In  Eq.  2.1-8,  the  3*n  output  matrix  B  selects  those  linear 
combinations  of  the  elements  of  x  that  are  being  estimated . 


The  following  discussion  provides  an  example  of  how 
this  formulation  is  used  for  gradiometer  test  data  reduction. 
For  convenience,  the  quantities  to  be  estimated  are  separated 
from  those  being  measured  when  defining  the  elements  of  x.  (I 
many  estimation  problems,  the  quantities  to  be  estimated  are 
included  among  those  being  measured  because  it  is  often  appro¬ 
priate  to  treat  measured  and  estimated  quantities  in  a  unified 
manner.)  For  GGSS  testing,  the  residual  gravity  disturbance 
(d(r)  =  6(r)  -  ^(r)]  is  to  be  estimated.  Therefore,  the  two 
conventional  gravity  quantities  needed  to  define  d(r)  [the 
point  disturbance  6(r)  and  its  local  mean  ^(r)]  are  listed 
first  in  x.  The  remaining  elements  of  x  are  used  to  represent 
the  actual  average  values  of  the  measured  quantities.  For 
gradiometer  data,  these  average  values  are  the  mean  gradients 
that  are  defined  by  the  appropriate  data  template  zones.  For 
example,  the  particular  prototype  template  used  in  this  report 
has  12  zones,  over  which  the  zz  component  of  the  gravity  gradi 
ent  is  averaged.  Therefore,  when  using  this  template,  the 
average  values  of  Tzz  for  these  zones  are  the  last  12  elements 
of  x  (these  mean  values  are  denoted  as  T  (i)  for  i  =  1  to  12) 


x  = 


6(r) 


Tzz(1> 

Tzz<2> 


T,_<12) 


(2.1-9) 


number  of  computations.  The  process  of  template  design  can  be 
performed  systematically  by  using  error  covariance  analyses. 
Therefore,  the  template  geometry  is  selected  to  be  not  only 
appropriate  for  estimating  the  residual  gravity  disturbance 
vector,  but  also  amenable  to  an  efficient  and  accurate  covari¬ 
ance  analysis. 

In  the  following  discussion,  the  measurement  vector 
z ^  is  a  list  of  mean  gravity  gradient  elements,  which  are  com¬ 
puted  during  the  data  processing.  To  define  these  quantities 
completely,  the  location  (including  altitude)  and  size  of  each 
averaging  window  must  be  specified. 

2.1.3  Error  Models 

The  last  part  of  formulating  the  estimation  problem 
is  to  model  the  errors  in  the  observed  data  and  the  a  priori 
uncertainty  in  the  values  of  the  gravity  disturbance  and  grav¬ 
ity  gradient  field.  The  observed  data  are  represented  by  the 
measurement  vector  of  appropriately  averaged  gradiometer 
measurements.  A  model  for  the  averaged  data  is  given  as  follows 

z2  =  H:x  +  y2  (2.1-7) 

In  Eq.  2.1-7,  v^  is  an  m* 1  vector  of  measurement  errors.  The 
n*l  truth  vector  x  represents  the  actual  values  of  the  gravity- 
related  quantities.  These  quantities  consist  of  those  that  are 
being  estimated  as  well  as  those  that  are  being  measured.  The 
m*n  observation  matrix  selects  from  x  those  linear  combina¬ 
tions  of  its  elements  that  are  being  measured . 

The  residual  gravity  disturbance  vector  d(r) ,  which 
is  to  be  estimated,  is  also  expressed  in  terms  of  the  truth 


vec  tor : 


POINT  IS  BELOW 
CENTER  POINT 


Figure  2.1-1  Data  Template  Windows  Superimposed 

on  Survey  Tracks  (Not  to  Scale) 

are  then  listed  in  a  vector  denoted  and  used  optimally  to 
estimate  the  residual  disturbance  vector.  The  locations  and 
sizes  of  these  subregions  are  selected  to  achieve  specified 
rms  accuracies,  while  minimizing  the  computational  effort . 

A  convenient  way  of  visualizing  this  processing  is  to 
imagine  a  template,  as  depicted  in  Fig.  2.1-1,  that  is  super¬ 
imposed  on  a  map  of  the  survey  tracks,  with  the  center  of  the 
template  located  over  the  position  where  the  gravity  departure 
is  to  be  estimated.  The  "windows"  in  this  template  define  the 
subregions  within  which  the  gradient  measurements  are  averaged 
The  shapes  and  locations  of  these  windows  are  selected  to 
achieve  specified  rms  accuracies  with  the  smallest  possible 


2.1.2  Measured  Gravity  Quantities 


During  testing  of  the  GGSS ,  the  measured  gravity  quan¬ 
tities  will  include  the  six  distinct  elements  of  the  gravity 

JL. 

gradient  tensor.  These  gradient  measurements  will  be  made 
along  nearly  straight  flight  paths  (within  ±100  m  of  the  nomi¬ 
nal  flight  path)  that  form  a  grid  over  the  survey  region.  For 
analysis,  all  of  the  measurements  in  the  survey  can  be  listed 
in  a  measurement  vector  z.  To  estimate  the  number  of  these 
measurements,  consider,  for  illustrative  purposes,^  flight  paths 
that  are  5  km  apart  and  form  a  square  grid  300  km  on  a  side. 

Such  a  survey  consists  of  122  paths,  each  300  km  long,  for  a 
total  track  length  of  36,600  km.  If  the  six  gradient  elements 
are  measured  at  1  km  intervals  (e.g.,  one  set  of  measurements 
every  12  seconds  provided  by  a  plane  flying  300  km/hr),  then 
the  measurement  vector  z  contains  219,600  numbers.  To  process 
z  optimally  (to  estimate  the  residual  gravity  disturbance  vec¬ 
tor  d(r)  at  the  center  of  the  survey  region),  219,600  simul¬ 
taneous  equations  must  be  solved.  These  computations  are 
impractical  because  of  the  size  of  z.  By  using  a  data  averag¬ 
ing  technique  to  reduce  the  dimension  of  the  measurement  vec¬ 
tor  z,  a  practical  estimation  algorithm  is  developed. 

To  estimate  the  residual  gravity  disturbance  vector  d(r) , 
the  gradient  measurements  can  be  averaged  over  selected  sub- 
regions  of  the  survey  area.  The  resulting  average  measurements 


*lt  is  recognized  that  only  five  of  these  are  independent. 
However,  the  redundancy  of  the  sixth  should  be  used  to 
improve  overall  accuracy. 

tThe  actual  test  flight  pattern  may  entail  a  greater  extent  and 
variable  track  spacing  which  increases  with  distance  from 
the  center. 
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In  Eq.  2.1-2,  the  local-mean  disturbance  vector  is  defined  as 
the  following  average  centered  on  position  r: 


where 


6  ( r ) 
— m  — 


ffv 


(&)  6 ( r-£ )  dp, dp. 


(2.1-3) 


(2.1-4) 


In  Eq.  2.1-3,  w(£)  is  the  weighting  function  at  shift  p. 
The  weighting  function  for  a  uniform  local  average  over  a  square 
geographic  region  centered  on  position  r  is  given  by  the  fol¬ 
lowing  equations  in  which  L  is  the  length  of  each  side  of  the 
averaging  area: 

w(p)  A  L'2,  -L/2  <  px  <  L/2 ,  -L/2  <  p2  <  L/2  (2.1-5) 
w(p)  4  0,  otherwise  (2.1-6) 


The  weighting  function  defined  by  Eqs.  2.1-5  and  2.1-6 
is  recommended  for  use  in  testing  because  (1)  it  provides  the 
needed  control  over  the  long-wavelength  variance  of  the  resid¬ 
ual  disturbance  vector,  and  (2)  it  yields  lagged  covariance 
functions  for  the  residual  disturbance  that  can  be  evaluated 
in  closed  form.  The  length  parameter  L  controls  the  long- 
wavelength  variance  of  the  residual  disturbance  d(r) .  The 
variance  at  wavelengths  longer  than  L  is  significantly  attenu¬ 
ated,  while  the  variance  at  wavelengths  shorter  than  L  is  only 
slightly  affected.  Because  the  lagged  covariances  of  d(r)  can 
be  evaluated  in  closed  form,  the  weighting  function  defined  by 
Eqs.  2.1-5  and  2.1-6  is  especially  suited  for  covariance  anal¬ 
yses  of  survey  accuracy.  On  the  basis  of  such  covariance 
studies,  algorithms  are  being  developed  by  TASC  for  processing 
the  survey  data  to  meet  specified  rms  accuracies. 


2-3 


A  prototype  template  for  estimating  the  z  residual 
disturbance  component  from  airborne  measurements  of  the  zz 
gradient  element  is  defined  by  the  14  x  8  documentation  matrix 
presented  in  Table  2.4-1.  This  template  is  a  prototype  example 
of  the  templates  that  will  be  used  for  processing  the  full  set 
of  gradiometer  test  data.  Later  templates  will  contain  more 
circuits  and  be  applied  to  additional  gravity  gradients.  The 
purpose  of  the  prototype  template  is  to  provide  an  interim 
example  of  the  type  of  data  averaging  and  data  documentation 
that  is  recommended  for  GGSS  testing,  as  well  as  to  offer  a 
rudimentary  data  processing  capability. 

TABLE  2.4-1 

DOCUMENTATION  MATRIX  FOR  PROTOTYPE  TEMPLATE 


ROW 

NO. 

ID 

xc 

(km) 

yc 

(km) 

w 

(km) 

A1 

(km) 

A2 

(km) 

A3 

(km) 

A4 

(km) 

1) 

1 

— 

0 

0 

0 

0 

0 

0 

0 

2) 

2 

-125 

125 

250 

0 

o 

0 

0 

3) 

7 

0 

0 

0 

0.6 

0 

6 

0. 

6 

0.6 

4) 

8 

-5 

5 

1 

0.6 

0 

6 

0. 

6 

0.6 

5) 

8 

-10 

10 

1 

0.6 

0 

6 

0. 

6 

0.6 

6) 

8 

-15 

15 

1 

0.6 

0 

6 

0. 

6 

0.6 

7) 

8 

-20 

20 

1 

0.6 

0 

6 

0. 

6 

0.6 

8) 

8 

-30 

30 

1 

0.6 

0 

6 

0. 

6 

0.6 

9) 

8 

-40 

40 

1 

0.6 

0 

6 

0. 

6 

0.6 

10) 

8 

-60 

60 

1 

0.6 

0 

6 

0. 

6 

0.6 

ID 

8 

-80 

80 

1 

0.6 

0 

6 

0. 

6 

0.6 

12) 

8 

-120 

120 

1 

0.6 

0 

6 

0. 

6 

0.6 

13) 

8 

-160 

160 

1 

0.6 

0 

6 

0. 

6 

0.6 

14) 

8 

-320 

320 

1 

0.6 

0 

6 

0. 

6 

0.6 

The  first  two  elements  of  the  truth  vector  x  described 
by  Table  2.4-1  are  the  point  z  gravity  disturbance  and  its 
mean  value  averaged  over  a  250-km  x  250-km  region.  Therefore, 
the  z  component  of  the  residual  disturbance  is  computed  from  x 
as  follows: 
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(2.4-1) 


dz(r)  =  Bx 

B  =  [1  -1  0lxl2)  (2.4-2) 


2.5  COVARIANCE  MATRICES 

As  explained  in  Sections  2.1  and  2.2,  to  implement  an 
optimal  estimator  for  a  particular  data  template,  or  to  com¬ 
pute  the  rms  accuracy  of  such  an  estimator,  two  covariance 
matrices  are  needed.  These  are  (1)  the  error  covariance  ma¬ 
trix  Pq  of  the  initial  estimate  Xq  of  the  truth  vector  x 
(Eqs.  2.1-15  and  2.1-16),  and  (2)  the  error  covariance  matrix 
Rj  of  the  errors  v^  in  the  measurement  vector  z^  (Eqs.  2.1-7 
and  2.1-18).  In  this  section,  these  covariance  matrices  are 
defined  for  the  prototype  template  defined  in  Table  2.4-1. 
Numerical  values  for  these  covariances  are  computed  using  the 
flat-earth  version  of  the  AWN  worldwide  gravity  disturbance 
model  (Ref.  2)  and  a  white-noise  error  model  for  the  gradi- 
ometer  measurement  errors. 

2.5.1  Gravity  Error  Covariance  Matrix 

With  reference  to  Table  2.4-1,  the  first  gravity  quan¬ 
tity  has  ID  No.  =  1,  and  is  the  point  value  of  the  z  gravity 
disturbance.  This  is  the  vertical  gravity  disturbance  at  the 
point  with  coordinates  x  =  0  and  y  =  0.  Therefore, 

x(l)  4  6z(0)  (altitude  =  0  km)  (2.5-1) 

According  to  Table  2.4-1,  the  second  gravity  quantity 

has  ID  No.  =  2.  The  table  indicates  that  x(2)  is  the  mean 

value  of  6  (r) ,  averaged  over  a  250  km  x  250  km  square  area 
z  T 

with  center  at  x  =  0  and  y  =  0.  Therefore,  with  r  =  [x  yl  , 
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(altitude  =  0  km) 

(2.5-2) 

The  third  quantity  documented  in  Table  2.4-1  has  ID 
No.  =  7,  which  denotes  the  point  value  of  the  zz  disturbance 
gradient  Tzz •  The  point  value  is  at  location  x  =  0  and  y  =  0 
and  an  altitude  of  0.6  km.  Therefore, 

x(3)  4  Tzz(0)  (altitude  =  0.6  km)  (2.5-3) 

The  remaining  gravity  quantities  documented  in 
Table  2.4-1  have  ID  No.  =  8,  which  denotes  mean  zz  disturbance 
gradients,  with  their  averages  taken  over  four  rectangular 
subregions  that  are  uniquely  identified  from  the  upper-left 
corner  coordinates  (xc,yc),  the  width  w,  and  the  altitudes  , 
A2 ,  A^  and  A^.  For  example,  the  fourth  quantity  is  defined 
as  follows: 


125 


x(2)  4 


(250)' 


JJ  6z(r)  dx  dy 


-125 


x(4)  4  (Ij_  +  I2  +  i3  +  I4)/4 


(2.5-4) 


5  5.5 

Ij  4  1/10  J  Tzz(£)  dx  dy  >  (altitude  =  0.6  km) 

-54.5  I 

(2.5-5) 


-4.5  5 

I2  4  1/10  !  /  /  Tzz(r)  dx  dyl  (altitude  =  0.6  km) 

5.5  -5  I 

(2.5-6) 
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i3  4  1/10 


(altitude  =  0.6  km) 


5  -4.5  ) 

/  /  TZz<£)  dx  d4 

-5  -5.5  ] 


(2.5-7) 


5.5  5 


I4  4  1/10  //  Tzz(r)  dx  dy  >  (altitude  =  0.6  km) 


4.5  -5 


(2.5-8) 


The  error  covariance  matrix  Pq  (for  an  initial  esti¬ 
mate  Xq  =  0)  has  been  computed  for  the  AWN  gravity  model  using 
the  above  definitions  and  their  extensions  to  the  other  grav¬ 
ity  quantities  documented  in  Table  2.4-1.  The  resulting  co- 
variance  matrix  has  been  provided  in  Ref.  1. 

2.5.2  Measurement  Error  Covariance  Matrix 

Previous  analyses  of  gravity  gradiometer  test  data 
(e.g.,  Refs.  4  and  5)  provide  the  basis  for  the  gradiometer 
measurement  error  model  used  in  this  report.  The  noise-like 
measurement  errors  are  typically  modeled  by  a  zero-mean,  white- 
noise  signal  component  added  to  a  random-walk  signal  component 
that  models  instrument  drift.  During  the  first  stage  of  proc¬ 
essing  the  raw  gradiometer  survey  data,  track-crossing  adjust¬ 
ments  will  be  applied  to  the  measured  gradients  to  suppress 
the  long-wavelength  errors  caused  by  random  instrument  drift. 
The  residual  errors  in  the  adjusted  data  will  then  consist 
primarily  of  random  errors  that  are  accurately  modeled  as 
zero-mean  white  noise. 


The  white  noise  in  the  individual  gradient  measure¬ 
ments  will  not  be  affected  by  the  track-crossing  adjustments. 

2  2 

Therefore,  the  mean-square  error  o  (E  )  caused  by  white  noise 
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in  a  single  gradient  measurement  can  be  computed  from  the  spec- 

2 

tral  density  S  (E  /Hz)  (two-sided  spectrum)  of  the  noise  and 
the  effective  averaging  time  t  (sec)  of  the  gradiometer 

aVc 

output  filter: 


2 

a  =  S/t. 


(2.5-9) 


For  example,  if  S  =  80  E  /Hz  and  tave  =  12  sec  (one  meas 

urement  every  kilometer  at  a  speed  at  300  km/hr) ,  then  the  rms 

error  of  the  individual  measurements  is  a  =  2.6  E  =  0.26  mgal/km. 

When  these  measurements  are  averaged  using  the  data  template, 

2 

the  mean-square  error  am  of  the  resulting  mean  gradient  is  de¬ 
termined  by  the  number  N  of  measurements  in  the  average: 


o2  =  a2/N 
m  ' 


(2.5-10) 


In  Eq.  2.5-10,  the  number  of  measurements  N  is  determined  by 
the  length  L  (km)  of  the  averaging  zone  in  the  data  template, 
the  speed  V  (km/hr)  of  the  survey  aircraft,  and  the  time  t 

aV  C 

(sec)  between  consecutive  measurements: 


N  =  L/(Vtave) 


(2.5-11) 


It  is  recommended  that  the  organization  of  the  gravity- 
related  quantities  in  the  truth  vector  x,  which  is  defined  by 
the  documentation  matrix  D,  should  also  be  used  to  organize 
the  measurement  vector  z^ .  To  achieve  the  same  organization 
in  both  vectors,  the  mxn  matrix  H^  in  Eq.  2.1-7  is  defined  for 
the  template  under  discussion  as  follows: 


H1  -  ^°mxp  Xmxm^ 


p  =  n  -  m 


(2.5-12) 


In  Eq.  2.5-12,  0  v  is  the  mxp  zero  matrix  and  I  is  the  mxm 

mxp  r  mxm 

identity  matrix.  This  definition  of  H^ ,  with  m  =  12  and  p  =  2, 
is  appropriate  for  the  example  template  defined  by  Table  2.9-1, 


in  which  the  first  two  gravity  quantities  are  used  to  character¬ 
ize  the  residual  gravity  disturbance  component  which  is  being 
estimated .  The  remaining  m  quantities  listed  in  Table  2.4-1 
and  the  truth  vector  x  are  the  measured  quantities  in  z^ .  The 
order  of  listing  these  quantities  is  preserved  because  of  the 
mxm  identity  matrix  in  Eq.  2.5-12. 

The  error  covariance  matrix  for  the  prototype  tem¬ 
plate  documented  in  Table  2.4-1  has  been  computed  using  a 

gradiometer  white-noise  model  with  two-sided  spectral  density 
2 

S  =  80  E  /Hz,  an  rms  error  per  gradient  measurement  of  a  =  0.26 

mgal/km,  an  aircraft  speed  of  V  =  300  km/hr,  and  an  averaging 

time  of  tave  =  12  seconds  per  measurement.  The  covariance 

matrix  is  diagonal  because  the  gradiometer  measurement  errors 

are  modeled  as  uncorrelated  from  measurement  to  measurement. 

2 

The  12  diagonal  elements  of  ,  expressed  in  (mgal/km)  ,  are 
listed  below: 


R1(l,l)  =  0.262  e  6.76  x  10'1 
R1(2,2)  =  0 . 262/40  e  1.69  x  10'3 
R1(3,3)  =  0 . 262/80  e  8.45  x  10-4 
R1(4,4)  =  0 . 262/120  e  5.63  x  10'4 
R1(5,5)  =  0 . 262/l60  e  4.23  x  10'4 
Rx(6,6)  =  0 . 262/240  s  2.82  x  10'4 
R1(7,7)  =  0 . 262/320  e  2.11  x  10'4 
R1(8,8)  =  0 . 262/480  e  1.41  x  10'4 
R1(9,9)  =  0 . 262/640  e  1.06  x  10"4 
R1(10>10)  =  0 . 262/960  e  7.04  x  10"5 
R1(ll,ll)  =  0 . 262/1280  e  5.28  x  10"5 


R1(12,12)  =  0 . 262/2560  e  2.64  x  10‘5 
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2.6 


ESTIMATION  AND  ERROR  COVARIANCE  EQUATIONS 


The  purpose  of  this  section  is  to  summarize  the  cal¬ 
culation  for  estimating  the  z  component  of  the  residual  dis¬ 
turbance  vector  using  Eqs.  2.2-1  to  2.2-7  and  the  data  template 
and  covariances  defined  in  Sections  2.4  and  2.5. 

For  an  initial  estimate  Xq  =  0  of  the  truth  vector, 
the  optimal  estimate  of  the  residual  vertical  disturbance  dz(r) 
(located  at  the  point  beneath  the  center  of  the  data  template) 
is  given  as  follows: 


d  (r)  = 
z  — 

G  z. 

(2.6-1) 

G  = 

DO 

p* 

t— * 

(2.6-2) 

B  = 

11  -1  olxl2) 

(2.6-3) 

K1  = 

P0HT'H1P0HT  +  Ri'_1 

(2.6-4) 

Hl  = 

l012x2  X12xl2] 

(2.6-5) 

In  Eq.  2.6-1,  z^  is  the  12x1  vector  of  averaged  gradiometer 
data,  averaged  using  the  data  template  documented  in  Table  2.4-1 
The  measurement  error  covariance  matrix  R^  is  given  in  Sec¬ 
tion  2.5.2  and  the  initial  gravity  covariance  matrix  Pq  has 
been  provided  (Ref.  1). 

The  error  variance  of  the  estimate  in  Eq.  2.6-1  is 

Pdd  =  El<dz(r)  -  Vi))2!  =  BP1bT  (2.6-6) 

Pi  =  Pn  - 


k,h,pa 


(2.6-7) 


1  /2 

For  the  numerical  values  provided,  the  rms  error  (P^  )  of  the 
estimate  of  the  residual  vertical  disturbance  is  3.6  ragal,  as 
computed  using  Eq.  2.6-6.  Other  data  templates  are  being  devel¬ 
oped  at  TASC  to  meet  the  accuracy  requirements  for  processing 
of  the  GGSS  test  data.  These  templates  will  be  described  in 
future  reports. 


3. 


SUMMARY  AND  CONCLUSIONS 


3 . 1  SUMMARY 

This  report  has  presented  a  practical  methodology  for 
processing  airborne  gravity  gradiometer  data  to  estimate  gravity 
disturbance  vectors  at  the  surface  of  the  earth.  The  methodology 
consists  of  two  stages.  During  the  first  stage,  the  gradiometer 
measurements  are  averaged  over  subregions  of  the  surveyed  area 
using  data  templates.  The  purpose  of  this  first  stage  is  to 
reduce  the  dimensionality  of  the  computations.  During  the 
second  stage,  the  averaged  gradiometer  measurements  are 
weighted  optimally  and  summed  to  estimate  the  gravity  disturb¬ 
ance  vector  (or  in  the  case  of  GGSS  testing,  to  estimate  the 
residual  gravity  disturbance  vector).  The  data  templates  are 
selected  to  meet  specified  rms  accuracy  requirements,  while 
avoiding  unnecessary  computations. 

As  an  interim  example  of  the  type  of  data  processing 
and  documentation  that  are  required  for  handling  real  survey 
data,  a  prototype  data  template  is  defined.  Associated  with 
this  template  are  gravity  and  measurement  error  covariance 
matrices,  which  are  needed  to  design  optimal  estimators  for 
this  template  and  analyze  their  rms  errors. 

Special  attention  is  given  to  the  unique  requirements 
for  testing  of  the  GGSS.  Because  the  limited  geographic  extent 
of  the  survey  area  makes  long-wavelength  accuracy  verification 
imprecise,  a  specially  defined  gravity  disturbance  is  intro¬ 
duced,  called  the  residual  gravity  disturbance  vector  (the 
departure  of  the  disturbance  vector  from  its  local  mean  value). 


The  residual  gravity  disturbance  vector  should  be  estimated 
during  testing  of  the  GGSS ,  rather  than  the  gravity  disturbance 
vector,  because  the  long-wavelength  content  of  the  residual 
disturbance  vector  can  be  more  easily  controlled. 


3.2  CONCLUSIONS 


Based  on  the  results  of  this  study,  the  following 
principal  conclusions  are  reached. 


•  A  practical  methodology  is  available  for 
processing  airborne  gradiometer  data  to 
estimate  gravity  disturbance  vectors  (or 
residual  gravity  disturbance  vectors  for 
testing  of  the  GGSS)  at  the  surface  of 
the  earth 

•  Data  processing  templates  can  be  de¬ 
signed,  using  error  covariance  analy¬ 
ses,  to  reduce  the  dimensionality  of 
the  data  processing  while  satisfying 
accuracy  requirements 

•  The  technical  approach  (minimum-variance 
estimation  based  on  averaged  measure¬ 
ments)  is  flexible  and  can  handle  mix¬ 
tures  of  gravity  quantities  appropriate 
to  future  survey  needs. 
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